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Abstract: We study anomalous charged fluid in 2n-dimensions (n > 2) up to sub-leading 
derivative order. Only the effect of gauge anomaly is important at this order. Using the 
Euclidean partition function formalism, we find the constraints on different sub-leading order 
transport coefficients appearing in parity-even and odd sectors of the fluid. We introduce a 
new mechanism to count different fluid data at arbitrary derivative order. We show that only 
the knowledge of independent scalar-data is sufficient to find the constraints. In appendix we 
further extend this analysis to obtain fluid data at sub-sub-leading order (where both gauge 
and gravitational anomaly contribute) for parity-odd fluid. 
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1 Introduction and Summary 

In past few years there has been much interest and progress in further understanding of rela¬ 
tivistic, charged, dissipative fluid in presence of some global anomalies. Presence of quantum 
anomalies play a crucial role in transport properties of fluid. The first evidence of quantum 
anomaly in fluid transport was holographically observed in [1, 2], The authors found a new 
parity-odd term (and hence a new transport coefficient) in the charge flavour current. The 
origin of this new term can be traced back to gauge Chern-Simons term in the dual super¬ 
gravity theory. Soon after these results were published, it was shown that the new parity-odd 
term in the charge current is essential because of the triangle flavour anomalies and the second 
law of thermodynamics [3] . In general the second law of thermodynamics (or equivalently the 
positivity of divergence of entropy current) imposes constraints on different transport coeffi¬ 
cients. The same constraint can also be obtained from the equilibrium partition function of 
fluid [4, 5]. Equilibrium partition function provides an alternate and a microscopically more 
transparent way to derive the constraints on these transport coefficients. A generalization of 
this approach for charged 17(1) anomalous fluid in arbitrary even dimensions up to leading 
order has been considered in [6]. 

In [7] Bhattacharyya et.al. studied parity odd transport for a four dimensional non-conformal 
charged fluid at second order in derivative expansion. In four spacetime dimensions the effect 
of anomaly appears at one derivative order and the parity-odd transport coefficients at this 
order are determined in terms of anomaly coefficient. In this paper the authors studied the 
transport properties at second order and found that out of 27 transport coefficients 7 are 
fixed in terms of anomaly and lower order transport coefficients. The goal of our current 
paper is to generalize this work to arbitrary even dimensions. In 2 n spacetime dimensions 
the leading effect of anomaly appears at (n — 1) derivative order. Hence the subleading 
corrections appear at nth derivative order. The aim of this paper is to study the constraints 
on transport coefficients appearing at subleading order. We innovate a systematic mechanism 
to compute different fluid data at arbitrary derivative order (parity odd or even). We list all 
possible scalars, vectors and tensors at any arbitrary derivative order in this paper. It seems 
to be rather difficult to find the independent sets. However, we argue that it is possible to 
get the correct constraint relations between transport coefficients even without knowing the 
independent sets of fluid data. 

Our analysis is not valid in two spacetime dimensions. In two dimensions the parity odd 
terms appear at zero derivative order itself, and hence parity-odd and parity even sectors are 
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not independent at any arbitrary order. Independence of these two sectors is important in 
our computation. 

In the parity-even sector, the leading correction appears at first order in derivative expansion, 
e.g. shear viscosity and bulk viscosity terms in energy momentum tensor etc. In this paper 
we have extended our calculation to include the sub-leading order correction ( i.e . second 
order corrections) to parity-even sector in constitutive relations in arbitrary even dimensions 
in presence of 17(1) gauge anomaly. This completes the description of fluid dynamics up to 
sub-leading order in derivative expansion (both in parity-odd and even sectors) in arbitrary 
even dimensions with abelian gauge anomaly. 

The organization of our paper is as following. In § 2 we explain our notation and perturbation 
scheme which we use in this paper. In § 3 we construct the partition function for both gauge 
invariant and non-invariant sectors and compute the constitutive relations from the partition 
function. We also describe the construction of the anomalous entropy current. § 4 is the 
most important section of this paper. Here we first describe how to construct fluid data at 
arbitrary derivative order. Next, we list all the leading and sub-leading order scalars, vectors 
and tensors which may appear in constitutive relations up to sub-leading order in derivative 
expansion both in parity-even and odd sectors. Although, we have not been able to find the 
‘independent’ parity-odd vectors and tensors at sub-leading order, this does not inhibit us 
from finding the constraints on the transport coefficients. We elaborate this issue in § 4.4. 
Finally, in § 5 we list the constraint on the transport coefficients up to sub-leading order. In 
appendices we explain the Kaluza-Klien decomposition (appendix (A)) and sub-sub-leading 
order counting (appendix (B)). 

2 Scheme and the Perturbative Expansion 

We consider a 2ra-dimensional spacetime manifold M^n) with metric ds 2 = G^dx^dx 1 ' and 
gauge field 1-form A = A^dx 11 . We want to study fluid dynamics in this background. A fluid 
is a statistical system in local thermodynamic equilibrium, which is generally characterized in 
terms of (covariant) energy-momentum tensor 'T p ‘ v , (covariant) charge current J^ and their 
constitutive equations 


= r p j P + ‘S', = x 


( 2 . 1 ) 


T = d*4. is field strength for A. Here we have introduced a U (1) anomaly 3 and a gravitational 
anomaly S '. The form of these anomalies is well known in literature [8]. Most of our work 
here will be concentrated on fluid upto subleading derivative order, where only 17(1) anomalies 
contribute: 


J = (n + l)<7 (2n) * T hn = (n + l)C (2n) 

~ 2 n 



( 2 . 2 ) 
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only starts getting values at subsubleading derivative order. Let us explain our notation 
here. 


• All the fluid quantities (like currents, transport coefficients, independent terms etc.) 
appearing in parity-odd sector, are denoted by ‘tilde’ (e.g. A). On the other hand we 
use no special notation for parity-even sector (e.g. A). Wherever applicable, A = A + A 
denotes the total quantity (parity-odd and parity-even). 

• V and V denote the covariant derivative and on M( 2 n ) and the equilibrium manifold 

respectively. We use A and * as wedge product and Hodge Dual on all manifolds, 
as no confusion is possible. 

Due to dissipative nature of fluid, it is not possible to write an exact generating functional W 
(or action) for fluids from which one can derive the energy-momentum tensor and charge 
current J Therefore we write their most generic forms, allowed by symmetries, in terms 
of fundamental fluid variables and their derivatives in a particular thermodynamic ensemble. 
In our analysis we consider the fluid variables to be temperature id, chemical potential 1 v and 
fluid four-velocity u M with = — 1. 

We prefer to work in Landau Frame , where all the dissipation terms are transverse to the 
direction of the fluid flow. Hence, we can decompose 'T fiv and J 11 as 

T 1 ™ = E{p,v)u*u v + fP 1 ', J 1 * = Q{d,v)u fi + T' 4 , (2.3) 

where IT 41 ' and are the most generic symmetric tensor and vector made out of fluid 
variables. In the Landau frame 

u,Ii> w = 0, = 0. (2.4) 

The easiest way to implement this is to project all vectors or tensors appearing in or LP 1 ', 
transverse to using the projection operator 

pru = QU" + u n u V' 


Since fluid is a low energy fluctuation about the local thermodynamic equilibrium, and 
T^ 1 can be expanded in derivatives of fundamental fluid variables 


pit 11 ' — tt^ 

11 _ li (0) 


+ n<r, + 


'Y'/x_ 


+•<<.) 


+’ f < 2 ) 


(2.5) 


where n| l yj and Y(( y . involves N number of derivatives on fluid variables. The terms on RHS 
can have the most generic form as, 


n (N) = E *)T Z )t + P^ E °(N)t (#, ^)S (JV)t , 


l w~7Tp c 




( 2 . 6 ) 


1 Actually v = p./'h, where y is the chemical potential. 
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where S^ N y, ^(N)t anc ^ ^-(N)t are a c °U ec tion °f all possible gauge invariant scalars, vectors 
and symmetric traceless tensors (collectively known as data ) respectively, made out of fluid 
variables and source fields at N derivative order. Yht corresponds to sum over independent 
terms at any particular derivative order. The data which is required for our computation has 
been enlisted in § 4. 

In eqn. (2.6), the expression for Il|((( r . and are fixed up to some undetermined coefficients 
appearing at each derivative order. Therefore, a fluid is characterized by an infinite set of 
such unknown functions ('T(jv)t, 0"(jv)t> jv)t), known as transport coefficients. Fluid up to a 
particular derivative order is characterized by a finite number of such transport coefficients. 
In general, these transport coefficients are not all independent. The second law of thermody¬ 
namics (or equivalently, positivity of local entropy current) imposes restrictions on different 
transport coefficients 2 [9]. Such relations among various transport coefficients are known as 
constraints. 

[4] uses a different mechanism to find ‘some’ of these constraints. The idea is to write an 
equilibrium partition function for the fluid and derive the energy-momentum tensor and 
charge current from that partition function. Because of dissipation it is not possible to write 
a generating functional (W) for the fluid. However, one can still write a generating functional 
in equilibrium configuration, which we denote by W eqb . Using W eqb one can find all the 
constraint relations involving transport coefficients which comes with data that survives at 
equilibrium. 

More precisely, if the theory has a timelike Killing vector we can write an Euclidean 
generating functional using the background fields and Killing equation on the decomposed 
manifold S' 1 x Here S 1 is the euclidean time circle along with time period and 

M(d— i) is the spacetime transverse to [4] has conveniently chosen = do- Therefore, 
one can decompose the background in Kaluza-Klein form, 

ds 2 = G fiU dx ti dx l/ = —e 2cr (d t + a^dx*) 2 + gtjdx l dx\ 

A = A(dt + affix 1 ) + Affix 1 . (2-7) 

For more details please refer appendix (A). Using this choice along with the Landau Gauge 
conditions and velocity normalization, the most-generic energy-momentum tensor and charge 
current in eqn. (2.3) on M(d) can be decomposed into scalars, vectors and tensors on S 1 x 

2 Similar restrictions are also applicable to non-relativistic fluids and has recently been addressed for charged 
non-relativistic fluids in [10]. 
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M(d- 1 ) : 


T ij = E(&, u)v i v j + ir ij , 

V = — e a \ E{d, v)v l \/\ + ViV 1 H- , V ^ 

\ v:i + 

f = e 2 " (^ ")(! ++ (ffS))' 

r = Q(#,v)v i + ?, 

J = — e a ( Q(i?, ^)y / l + ViV 1 H- . ^ ) , (2.8) 

V Vi + vi«v 

where 

T = Too, V = T l 0 , V 3 = V j ; J = Jo, T = T, 

and 

vf = n 00 , Tf* = n* 0 , TfV = TV 3 , = To, T = T®, u = no,^ = u\ 

Indices on .Ad^-i) are raised and lowered using g*- 7 . Details of Kaluza-Klein decomposition 
of fluid variables and background fields have been given in appendix (A). 

Since the fluid we are considering is in local thermodynamic equilibrium, we can write the 
fluid variables as a spatial derivative expansion about their equilibrium values 

t? = $ 0 + A^tf + + ... 

v = v 0 + A (1) n + A+ ... 

j = j 0 + A (1 V + A (2 V + .... (2.9) 

The terms with subscript ‘o’ are the equilibrium values, while A^-* designates the N th deriva¬ 
tive corrections 3 . The zeroth component of fluid velocity uq = v also gets derivative correc¬ 
tions which are determined by the corrections to v % using the four-velocity normalization. 
Similarly all the transport coefficients can also be expanded using the Taylor Series expan¬ 
sion 

a($, v) = a 0 {J 0 , v 0 ) + A^a + A+ .... (2-10) 

Therefore the energy-momentum tensor and charge current receive two fold derivative cor¬ 
rections. First of all we write these expressions as a derivative expansion in terms of fluid 
variables in eqn. (2.5). Secondly, each term in that expansion can be further expanded around 

3 In this paper, A ^ A denotes parity-odd n th derivative corrections to a fluid quantity A, while A 
represents the parity-even n th derivative corrections. Entire derivative correction is denoted by A^A = 
A ( ”>A+A W A 
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the equilibrium values of fluid variables according to eqn. (2.9). Thus we finally get 


? = 


tt 13 

n o( o) 
A 

?o(0) 


+ 

+ 


A(1) ^o) + Km 


+ 


+a < i, *; i 1) +< 3 (2) 


J + [A<%, + A<%, + f 0 


o(2) 


( 2 . 11 ) 


Expansion of time components can be determined from these using Landau gauge condition 
eqn. (2.4). 


We choose the equilibrium convention for $ and v by identifying their equilibrium values to 
be the red-shifted temperature and Wilson loop in the lower dimensional theory 


— — /3 0 — P\J—Gqq — (3e a , v Q — /T4o- (2-12) 


In the next section we construct the equilibrium partition function and obtain energy-momentum 
tensor and charge current in terms of background data following [4]. After that, we compare 
these stress tensor and current with the fluid stress tensor and current order by order in 
derivative expansion to find the constraints among the transport coefficients at any partic¬ 
ular derivative order. A typical constraint will connect transport coefficients at equilibrium 
{a 0 ($ 0 , z^,)} and their derivatives with respect to i9 0 and u 0 (up to a particular derivative 
order) 

f({ao(^To)},{9ao(^To)}) =0. (2.13) 

We can extrapolate this constraint to non-equilibrium configurations: 

^( {«($, v)} , {da($, v)} ) = 0, (2.14) 

while doing this, we are making an error of at least one derivative order higher, which will 
be compensated at next derivative order computation. This is how we find the generic con¬ 
straints among fluid transport coefficients. Please note that while the equality constraints 
determined by this procedure are generic, the inequality constraints are not determined by 
this method. 


3 Equilibrium Partition Function 


The equilibrium partition function 4 W eqb of the theory can generally be disintegrated into 
two parts: 

w eqb = + W[f y (3.1) 


4 The partition function may be thought of as the Euclidean action for the fluid living on the background 
with coordinate time t compactified on a circle of length {3 
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The first one is the ‘conserved’ partition function which is gauge and diffeomorphism invariant, 
and generates conserved part of currents denoted by 7^, J{p\ • The other piece is not gauge- 
invariant and is referred to be ‘anomalous’ partition function. It generates anomalous piece 
of ‘consistent currents’ which will not be gauge-invariant in general. By defining a consistent 
subtraction scheme (Bardeen-Zumino currents), we can make these anomalous currents gauge 
invariant (see [11] for details) which we denote by 7Their value at equilibrium is 
fixed by anomaly, and upto subleading order is given by: 


'(A) ~ 


'tH 

J {A) 


n , 

_ 2C (2„) £ »+lC„ +1 tf>r +1 * (»„ A A X^~ m p <> 

m— 1 
n 

_ 2C (2n) -+ 1 C m+1 ^C+ 1 Zj m) < ) , 

m =1 
n 

- C (2n) Y, (n + 1) n C m # 0 v™ * (u 0 A A X^ n ~ m) y 

m= 1 
n 

m=l 


(3.2) 


(3.3) 


Here {X\, X 2 } are {—fldu, d.4 + fli/du} projected transverse to u and their equilibrium values 
upon KK reduction reduce to: {/ 1 , / 2 } = s $da, dA Hence we have 5 : 


J U) 


3 (T q2 T71+1 7* 


C'( 2n ) ^ n+1 C r _ 

ra=l 

n 

_ c ,( 2 n) ^ ( n + 1) n C m ti 0 ls™l\ 


n 

'o(ra) ’ 

o' 

II 

II 

■os 

(3.4) 

i 

o(m) ’ 

J(A) = 0. 

(3.5) 


m=l 


Let us now concentrate on Wpy It’s variation on background (2.7) will determine the con¬ 
served currents: 

1^, 


8W$ = J d 2n x VG 


And hence, 


8Wft 
~ 5G, U ’ 


-§ TfflG^ + JfoSA, 


SWfft 

sffd _ V-') 

J ^ ~ 5 A, ■ 


(3.6) 


(3.7) 


Kaluza-Klein decomposition of eqn. (3.7) gives, 


SW?£ 

r r l ) _ o„Q T; 

g . 


<sw$ 


P [C) + ^oJlc) ~ ^ 


r <o = <^ ». 




5 To get these and some further results we have to used the ideal order results v 0 — —e a ,v l 0 = 0, which we 
will derive in § 5. We use it here to simplify the notation. 










(3.8) 


J(C) ~ $°~ 


sw; 


eqb 


(C) 


SA; 


J iC) = -e° 


SW, 


eqb 

(C) 


SVn 


Here we have switched the basis to d 0 = e _<T //3 and v 0 = (3A for later convenience, a 1 is 
the Kaluza-Klein gauge field. Note that while Wf^ is gauge invariant, its integrand does not 
need to be. We can include a typical Chern-Simons term to it, which is defined such that its 
integral is gauge invariant 6 : 


r2n—1 


= - d 


' 2 ra —1 


' M(2n_ 


Xyfg 


(2n —1) 



Cm-iCm-iAil 0 ( m j + dC n a,il 0 ^ 


(3.9) 


Here C m ’s are constants. This is indeed a valid Chern-Simons form as at equilibrium is 
just made of Chern classes of /i and f 2 - 


fym) = * 


A(m—1) j.A (n—m) 
1 A J 2 


(3.10) 


For the gauge-invariant integrand, we assume that curvature scales of the background 
is much much larger than the mean free path of the fluid, therefore the whole manifold can 
be thought of as union of various flat patches. The system can be thought of in thermal 
equilibrium in each local patch. On each patch we can define the euclidean partition function 
locally, hence giving us 

W [c) = J d 2n ~ 1 x y /g (3(x)P(x) + J d 2n - l x^*I 2n ~\ (3.11) 

where P(x) is local thermodynamic pressure and /3(x) is local thermodynamic temperature. 
Given pressure, we can use the thermodynamic relations in local patch 

dP = e ^J-d r d + dqdv, e + P = ’ds + uq, (3-12) 

v 

to define energy density e, entropy density s and charge density q of the fluid. All are functions 
of d and v. We can expand W^ around its equilibrium value as 

= J d 2n - x xj~g p o Po + A W^ b y (3.13) 

Derivative correction to the ideal fluid partition function is denoted by A Wf^, which will 
contain all the possible gauge invariant scalars made out of background metric and gauge 
field components at a particular derivative order. We have computed these scalars (till the 
derivative level of our interest) in § 4. 


Collating together the conserved currents in eqn. (3.8) and the anomalous pieces in eqns. (3.4) 
and (3.5), and varying the Chern-Simons terms in W^ i.e. eqn. (3.9), we can finally 
write: 


T ij = 2iV 


sw: 


eqb 


(C) 


$9i 


6 We have left the terms in I 2n 1 


which can be related to others upto a total derivative. 
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f l + e a d 0 u 0 J l = do 


5W, 


eqb 

(fl) 


Scti 


A £ + »<?„*• w 

m =1 ' ' ' 


swf. 


r eqb 

'(C) 


,r = ■do^r^- - do n_1 c 


T — »2o- o2 (C) 

° <h? 0 ’ 

n(n + 1) 


m—1 


771=1 


m 




J =-e° 


5W, 


eqb 


Cc) 


Sv 0 


(3.14) 


Comparing these to the most generic fluid expressions in eqn. (2.8) we can compute the 
constraints. Thus, we see that it is only the gauge invariant part Wf^ of the partition 
function that we need to evaluate at any desired order. 


3.1 Anomalous Entropy Current 

In last section we reviewed a procedure to get equality type constraints among fluid transport 
coefficients. It is generally known that these very constraints can also be get by demanding 
existence of an entropy current whose divergence is positive semi-definite. The most generic 
Entropy Current can be written as: 

J l S = J'kc) + ^S(Ay (3.15) 

where ^ is the part which captures the explicit dependence on anomaly coefficients. How¬ 
ever, the other piece J$(q) can get implicit dependence on the anomaly coefficients through the 
fluid equations of motion. We need to demand this current to be positive semi-definite, 

^ = Vm^(o) + ^ ^s(A) — 0’ (3.16) 

whenever EOM are satisfied. For equilibrium fluid configuration, both the pieces can be 
demanded to be positive semi-definite separately. Such decoupling is not always possible, as 
the fluid equations of motion depend on anomaly coefficients, which can induce some implicit 
anomaly dependence in ■ However, for equilibrium fluid configurations, the equations of 
motion are trivially satisfied and thus entire information of anomaly can be incorporated in 
Jg( A) ■ Hence, if any part of V n<Jg( C ) cou Pl e 1° ^ ^s(A )> ^ ie res P ec tive transport coefficients 
will be determined in terms of anomaly coefficients, and hence will be present in 
first place. Therefore all the information about constraints among fluid transport coefficients 
is encoded in the existence of Jg^ C y In [12, 13] the author gives an explicit construction of 
entropy current from Eqb. Partition Function. 

Now concentrating on the second term: at equilibrium, V /JA 7^(/q > 0, since it does not have 
any independent coefficients, just constants, one cannot apply any constraints for it to be 
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satisfied. Therefore must be exact, 

exact terms, and hence we can choose J$( A ) 
hydrodynamic frame': 


But any current is always ambiguous upto some 
= 0 equally well. We can hence write in a generic 


Va) = 


-7 

J (A) 


(u v) 

(A) U (C) 


+ 2g, 




(A) U (C) 


icy 


+ 2 q'/^u 


*0 

{A) 


+ Kv 




<1 U '(A) + T(!A)> 


0 - SU 


(3.17) 

(3.18) 

(3.19) 


Note that in the expression for 7^ we have used the fact that anomalies are parity-odd. Now 
depending on the choice of hydrodynamic frame, these conditions can be used to determine 
anomalous dissipative parts of the various currents. For example, if we define such that it 
does not contribute to anomaly, i.e. = 0, we will get: 

<A) = -T(A) U ^Y Ka) = . %)=»• 

(3.20) 

This is the neatest frame for anomalies. Similar results for U( 1) anomaly were derived in 
[14], however these expressions are also applicable to gravitational anomalies * 8 . Here we 
present explicit expressions for the anomalous parts of currents, in presence of both U( 1) and 
gravitational anomaly. Following the generic expressions given in [8], these can be computed 
directly from the anomaly polynomial. The anomaly polynomial in 2 n dimensions up to 
(n + 1) derivative order is given as [15], 

V = C'( 2n )j rA ( n+1 ) + A Trpt A 91], (3.21) 

where, C^ 2n ^ is gauge anomaly coefficient which we have already introduced in the last section 
and c m is gravitational anomaly coefficient. The two form 91 is defined in terms of the Riemann 
tensor as, 

9 \ a p = TZ a ^ s dx' y A dx 5 . (3.22) 

Taking appropriate derivative of the above, one can find explicit expressions for anomalous 
parts of the currents. The leading part of the currents proportional to the gauge anomaly co¬ 
efficient C( 2n ) have already been given in eqns. (3.2) and (3.3). Here we present the subleading 
order contributions to currents coming due to the gravitational anomaly, 


— c 
J {A) ~ 1 

n— 2 


>(77,-1) *(u 0 AjF A ( n 2) AA aji 


- 29T /3 


+ ^2 n - 2 C m (y 0 V 0 ) m * (uo A C/ A(m - 1} A JTA(n-2-m) A ^a/3 _ A _ A ^[7))' 


m=l 


(3.23) 


'We have used the thermodynamic functions e, q,s here, which will be explicitly proved in § 5. 

8 restricted to equilibrium configurations 
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where, 


A fiu = 7 ; ( U fW - 4—u 0 [j,P M ] a V a t? 0 ) , U, w = 2P[ fia P u ] /3 'V a u^ 


(3.24) 


The heat current has the form, 
1 


I (A) — Cm 

n— 1 




★ a j' A ( n_i )y i A"^A a/3 


+ n- 1 C m (rn - l)(tf 0 Z/ 0 ) m * («o A C/ A ( m - 2 ) A jrA(n-l-m) A ^a/3 _ A _ A af} U 


m=2 


n— 1 


-2 ^ n ~ l CmiVoVoT * (u 0 A A JTA(n-l-m) A A afi _ A 


m= 1 

Finally, the stress tensor looks like, 

~n— 1 


. (3.25) 


T^ v — 4c V 

> — ^ L '-m v p 


Y n_1 C m (^ 0 z/ 0 ) m * (u 0 A P A ( m - 1 ) A A (w<r - hf^uyfi 

-*(u 0 AF A ('*- 1 )A p (j'‘ ) -2^'qJ^. (3.26) 


_m=l 


Instead if we are working in Landau Frame, where q ^ = 0, we will get condition: 


- T(Z )Uu{c) = (eG^ + T,Z) u v(A) = (eG^ + HZ 


(C )) U PA) 


(c)j u HAy 


(3.27) 


We need to invert (eG^ v + 11^ ^ , which can be done perturbatively in derivatives. To leading 
order: 

1 M " (3.28) 


U (A)~ e + p^{A) u PC) + • • 


and hence 


Th, = j ( % + Zf>T(% u W + - ■ ■ ’ fi f. A) = 2^7 (<7 + 7^Ao)“»(c) )+■■•. (3.29) 

(3.30) 


%) = ApPa^c) + • • • = It7“»(C) + - ■'Tfo + • • • • 


As showed by [14], in presence of just 17(1) anomaly, it gives the exact result of Son-Sorowka 
[3]. To write a similar expression for gravitational anomaly in Landau Frame, one will need 
to find anomalous velocity to subsubleading order, which might be non-trivial. 


- 12 - 










4 Counting of Independent Terms 


This section is dedicated to develop a systematic procedure to compute independent fluid 
data (vectors, tensors transverse to velocity and scalars). First we will review the counting in 
parity-even sector in generic dimensions. Then we will extend this idea to parity-odd sector 
in generic dimensions at arbitrary derivative order through a procedure we call ‘derivative 
counting 

After describing the generic procedure, we explicitly construct leading and sub-leading order 
parity-odd and even terms which are important for our current work. Many of these terms 
vanish in equilibrium. In tables tables (1), (2) and (5) we list all the leading and sub¬ 
leading terms both parity-odd and even and check if they survive at equilibrium. Further 
in appendix (B) we extend this counting procedure to parity-odd subsubleading derivative 
order fluid. For this reason we will keep our illustrations in the construction explicit to 
subsubleading order. 

4.1 Parity-even Counting 

In this subsection we present the parity-even counting in generic dimensions. One can always 
count independent data in the local rest frame (LRF) of the fluid, which turns out to be 
easier. We can later covariantize the terms to a generic reference frame by following simple 
(and generic) rules 9 . In LRF, the fundamental quantities are 

• Temperature - •&, Chemical Potential - v. 

• Derivatives of fluid velocity 10 - dou l , d^u 1 . 

• Field Tensor - T l ' J , £ i = T iv u v . 

• Curvature - VJ /jk{) . lZ l0k0 . 

All other quantities are merely derivatives of these fundamental quantities. Since LRF is 
locally flat, we are using the coordinate derivatives d Q and <9j. We introduce a notation for 
parity-even terms which will be useful later in parity-odd counting. Terms with d derivatives 
and i indices will be denoted collectively as (^ ,i,d ). When working at equilibrium, it is also 

9 The rules can be summarized as: replace 1) ‘0’ index with contraction with iA, 2) %’ indices with a 
projection along P ^ v , 3) d with V, 4) w ith e ni v u* 2 v 2 ---i*nvn u ^n anc j g na gy 5 ) p U t all extra factors 

of projectors and velocities on left-most, so no derivatives act on them. 

“«% = — 1 would imply u^d = 0 and hence in local rest frame duo = 0 . 
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convenient to define 11 : 


(4.1) 

(4.2) 


S pv = 2V^u v \ U’ jv = 2V^u v \ 

XfP = {-dlA"h pi ,a P l ' 6 X a 8 + 01/*/'"'} , A = 1, 2. 

The purpose of above notation is revealed in Kaluza Klein formalism: at equilibrium only 
spatial components of Xj^ survive which land exactly to fi defined by: 

a A = {#a\ A 1 } , ft = VV A - V j a\. (4.3) 


In the same spirit we define 


]C puprT = P pOL P v0 P^P a5 lZ OLfil5 - (u pu U p(J + ^U lip U 1/a - ^U up U' l A . (4.4) 


Only spatial components of survive at equilibrium, and they exactly match RP kl . The 

usage of index ‘A’ is purely to facilitate counting and computations. Similarly we define 

= {&■,»}■ 

Bianchi identity: 

In counting, we will extensively use the Bianchi identity to get rid of many terms, so it would 
be worth to spend some time on it. The Bianchi Identities for Field Tensor, Vorticity and 
Riemann Tensor take the form: 


^ [pP~ up\ V [pP'i/p]crS '^'[pup]a O' (4-5) 

However our redefined variables X\ and /C do not satisfy Bianchi Identities. But nevertheless 
we can always use these identities to relate 


^ [ll^Aup] J 


^ [pPvp]a51 


K 


[pvp\(Ti 


(4.6) 


to other terms, and hence we can safely get rid of these in the following computation. In rest 
frame especially (or at equilibrium in any generic frame), one can check that X\ and /C also 
satisfy Bianchi Identities. 


Killing equation: 


11 Our conventions are: 

A lp " ] = ^P p a P u p (A ap - A 0a ) , 

A lij] = ^ ( A ij - A ji ^ , 


A^ v) = (A aP + A Pa ) , A kpv) = A^ v) - J^P a pA' 

A w) = ^ ( A ij + A ji ^j , A {ij) = A (ii) - -gijA ij . 
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If the theory has a unit Killing direction we have the following Killing equation for a 
general tensor 


£ujT aia 2 - = q 


u^V„T aia2 - 


y ^ 'j'Oii 
k 


...otk-ivak+x-y aUJ ak , 


(4.7) 


which in local rest frame becomes 

d 0 T a fo- = 0. (4.8) 

Therefore if we are considering a theory at equilibrium, we do not have to consider the do 
derivatives. Secondly, the Killing equation for metric G is given by 

V^ Q + V°c/ = 0. (4.9) 

Taking “ 7=7 = and using Killing Equation for scalars this translates to: 

VV* + V“r/ = 0 . (4.10) 


Hence in local rest frame S ’*- 7 = d l u 7 + d^u 1 = 0. 


4.1.1 First Derivative Order 


Below, we compute all possible terms at first derivative order in LRF. 

1 . ( 2 , 2 , 1 ): 


2 . ( 1 , 1 , 1 ): &■& A , dou* , £ i 

3. ( 0 , 0 , 1 ): S\, d 0 Va 


However all these first derivative terms are not independent on-shell. Using first order equa¬ 
tions of motion one can eliminate some of them. The equations of motion are given by 
eqn. ( 2 . 1 ) (at equilibrium) 


1 . ( 1 , 1 , 1 ): d/T^ = F a j a + % u 

2 . ( 0 , 0 , 1 ): d^o = 9^ = 2 - 


Using these equations we have killed the boxed terms in the counting. 


4.1.2 Second Derivative Order 

Below we list all possible pure second derivative terms. By pure we mean they are not product 
of two first derivative terms. Product of two lower derivative terms are called composite 
terms. 

1. (2,4,2): K,d kl 
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Table 1: Independent Leading Order Parity-even Data 


Name 

LRF 

Covariant 

Equilibrium 

0 

1 oi 

2° i 

1 CM 

2° M 

0 

K 

<9*1? A 

P /jq V q 'i? a 

V*l?Ao 

v? 

S i - t?iR 2 * 

- tfiKf 

0 


\ S&) 

Ig^v) 

0 


2. (|,3,2): &S' jk , d'wf, U ijk Q 


3. (1,2,2): 


d 0 S ij 


doX? , do?%, TV^ 0 , IC m { 


4. (i, 1,2): &do0A, dodou* ,d 0 S i ,d i S i l,d i Xji,TV a l 


aO 


5. (0,0,2): 


d 0 s\ 


dodo# A , djd l d A , 5^*, K ab ab , ft a 0a0 


Here also all the terms are not independent because of equations of motion. The second order 
equations of motion are given by, 


1. (1,2,2): d&dpTrt = d^ k (F^ a J a + T*)), d^d^T^ = d^ k (F^ a J a + T*l) 

2. (|, 1,2): dod„r> n = d 0 {F a J a + % l ) , <9* <9^ = -& ( £ a J a - % a u a ) , 

d l d^J^ = <9*3 

3. (0,0,2): did^T* = di {F a J a + £>), = -d 0 (5‘ Q - T Q u Q ), 

dod^J 11 = < 9 0 3 


Again we have killed boxed terms in the counting using equations of motion. We have 
provided a list of all terms till second order (also composites) in covariant form and their 
equilibrium values in tables (1) and (2). We can iterate this procedure to further derivative 
orders as required by the cause. Note that, for a pure term at IVth derivative order, the 
maximum number of indices possible are iV + 2; we will need it later. 


4.2 Parity-odd Counting 


In this section we shall compute the parity-odd leading and sub-leading derivative fluid data. 
Calculation in parity-odd sector is a lot more cumbersome, even in LRF. We introduce here 
a scheme called ‘derivative counting’ to compute these terms step by step. Any parity-odd 
term in (2n)-dimension must have a (2 n — l)-dim Levi-Civita involved in LRF 




(4.11) 


We are interested in constructing all possible scalars, vectors and symmetric tensors using 
it. A bit of thinking will reveal that one needs at least (2 n — 2)-rank parity-even tensors to 
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Table 2: Independent Subleading Order Parity-even Data 


Name 

LRF 

Covariant 

Equilibrium 

SiA 


P af) V a VpV A 

V'V^Ao 

^2(Ar) 

d l d\di'&r 

r 

V^AoVitfro 

S 3 (Ar) 

xzx rij 


fZMj 

s 4 

tc 

/c 

R 

s 5 

diVi 


0 

s 6 

Koo + 

— — y -tfiX-p 1 Xuj 

uWKfu, + |S U 
—2^82,11 — 3 ;A® 3 .ii 

0 

S7A 

VkV'M 

RVV 

0 

s 8 

viv 3i 


0 

S9 

02 

©2 

0 

S10 

S i3 Sij 

S^S, W 

0 

v M 

v 1 A 

dkXj? 

P^P a pV a X^ 

V fc /f 

V^r 

XfV Yk 

xrvra 

/a V fc dro 


do Vi 

P^u a V a V 3p 

0 


diS i3 

P^\ 7 a S a p 

0 

v t 

n\ - ±d k x ki - J*. x[ k d k d 

P^u a K a p - - ^Vi, n 

0 

v M 

V 6A 

QVi 

erf 

0 


QVi 

evy 

0 

V^ A 

s ij v Aj 

S^Vav 

0 


S ij v 3j 


0 

v M 

v 10 A 

xgv 3j 

xrv*, 

0 

rr\f-Ll' 

J - 1 A 

d^d^dA 

P^ a P^V 0 Vpi 9 A 

V^V^Ao 

rp/XZ'' 

± 2(Ar) 

dH K d j ^v 

p{/x*pv)Py a 0 aV^ 

V^AoV^tfro 

rp/XZ'' 

± 3 (Ar) 

■y{ik -yj) 

X Yk 

y(H<% "V^) 

fWffi 

2 A JTA: 

rp/iZ' 

± 4 

/C<b> 

K^ w) 

R(H) 

rp/il'' 

± 5 



0 

rp/iZ' 

^Vo + 

p^p u )°u a u Pn pa(J p + 

0 

± 6 

- 2 ^a<w>0 - 

O 1 rp/ZZ' 1 rp/ZZ 2 

— Z ^ J - 2,11 40 * A 3,11 

rp/iZ' 

± 7 A 

y( l yd 

V A v 3 

Rdf 

0 

rp/iZZ' 

± 8 

rr 

Rdf 

0 

rp/il'' 

± 9 

©o-d 

©cR"' 

0 

rp/ZZ^ 

- L io 


S^ a S u) a 

0 

rp/ZZ' 

11 A 

s^Y,, 

S^ a X^ a 

0 


be combined with e n2 - ?2 '" ln - 7n for this purpose. One can subsequently form a list of parity-odd, 




types: 

V e : Vectors with free index on e (2n — 2 rank parity-even tensor contracted with e). 

S: Scalars with all indices contracted with e (2 n — 1 rank parity-even tensor contracted 
with e). 

T e : Tensors with one free index on e (2n — 1 rank parity-even tensor contracted with 
*)• 

V f. Vectors with free index not on e (2 n rank parity-even tensor contracted with e). 

Vp: Vectors formed of contraction of two non-e indices with free index on e (2n rank 
parity-even tensor contracted with e). 

T/-: Tensors with no free index on e (2 n + 1 rank parity-even tensor contracted with e). 

S*^: Scalars formed of contraction of Tj (2n +1 rank parity-even tensor contracted with 
*)• 

Tp: Tensors formed of contraction of two non-e indices with one free index on e (2 n + 1 
rank parity-even tensor contracted with e). 

V^: Vectors formed of contraction of two non-e indices with one free index not on e 
(2 n + 2 rank parity-even tensor contracted with e). 

\CC : v ec tors formed of contraction of four non-e indices with free index on e (2 n + 2 
rank parity-even tensor contracted with e) 


and so on. 

Here we note that given D derivatives, one cannot construct a parity-even term, pure or 
composite, with more than 2D indices, because (2,2,1) and (2,4,2) have the highest index 
to derivative ratio, which is 2. Therefore, if we are interested in a fluid at (n — 2 + s) 
derivative order (s = 1 corresponds to parity-odd leading order and so on), we can get at 
most 2(n — 2 + s) indices. The list of parity-odd data types we gave above is complete till 
subsubleading derivative order (s = 3). 

Independent Data Types 

We should emphasise that not all parity-odd data-types listed above are independent. The 
dependence comes from that fact that when we are working in 2n — 1 dimensions, any antisym- 
metrization over 2 n or more indices will vanish. Given that we are dealing with parity-even 
tensors of arbitrary rank which are to be contracted with e, there are a whole lot of these 


data 

1. 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 

10 . 
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antisymmetrizations possible. Hence, to find the independent data-types becomes highly 
non-trivial. 

Let’s look at a special case of this dependence. We construct a 2n-antisymmetrization, 


*l—*2n-l 


= 0 , 


(4.12) 


therefore, 

2n—l 

ii...i 2n -i Akik2—kt _ \ ' /i\o+l,feiii...ian-2 Axkz—kt 

*l-*2n-J Z_^ ' 1 *l...*o-l**o-*2n-2‘ 

a=l 


(4.13) 


The consequence of this is that the data types [ ]/ (i.e. ones with a free index not on e) can 
be expressed in terms of [ ]p (i.e. the ones with a free index on e and an extra contraction). 
Hence data-types [ ]/ for example V/,T/,Vj are not independent. 

Note that this result is only based on a specific form of 2n-antisynunetrization (eqn. (4.12)). 
One can in principle go on with any random antisymmetrizations over 2 n or more indices and 
find relations among the data, which as it turns out, is not a trivial task to do. We will come 
back to this issue in § 4.4. For now we continue with the counting. 

4.2.1 Derivative Counting 

We have classified parity-odd terms in data-types based on the number of parity-even indices 
required. We want to construct all allowed parity odd terms with D derivatives. We ob¬ 
serve that it is not required to include all parity-even data type of the form (r, i, d ) in this 
construction. We will show this below. 

For a parity-odd fluid at D = (n — 2 + s) derivative order, we need to construct all the D 
derivative parity-even terms with number of indices ranging from 2D (the maximum possible) 
to 2 (D + 1 — s) (= 2n — 2, the minimum required), i.e. 

2 (D + 1 — s) < No of indices of a parity-even D derivative term < 2D. 

These D-derivative parity-even terms can be constructed out of pure derivative terms. We 
need not consider pure terms with self contractions in parity-even data types as they have 
been included in our counting procedure. 

We now want to argue that not all parity-even data-types are required for this construction. 
For a data-type (jr,i,N) to be included at least once, the following combination with (2D — 
2N + i) indices must be included: 

(D-N) x (2,2,1 )®(±,i,N) 
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Table 3: Parity-even Data-types - Surviving at Equilibrium 


Data Type 

Decomposition 

Local Rest Frame 

Equilibrium 

(2,2,1) 



fj 

J A 

(1,1,1) 


# l?A 

V^Ao 

(2,4,2) 


fcijkl 

fiijkl 

(|,3,2) 

#(2,2,1) 

d l Xf 


(1,2,2) 

#(1,1,1) 

##0 A 

V ! V#ao 

(|,5,3) 

#(2,4,2) 

gij^jklm 

yiftjklm 

(|,4,3) 

##(2,2,1) 

##T“ 

V^'/a 

(1,3,3) 

##(1,1,1) 

###'<5 A 

V*VlV fc 7?Ao 

(|,6,4) 

##(2,4,2) 

gi gjj^klmn 

V*V J i# mrl 

(|, 5,4) 

###(2,2,1) 

###T( m 

v*v j v fc /| m 

(1,4,4) 

###(1,1,1) 

####$ A 

V*VlV fc V^ Ao 

(b 7,5) 

###(2,4,2) 

###/C Zmno 

\/ i \/J\7 k Jl lrnno 

(1,6,5) 

####(2,2,1) 

####<% 

V*vlv fc v'/^ n 

(|,8,6) 

####(2,4,2) 

# # # # fc mno P 

V*V J V fc V z R mnop 


Since the minimum rank of this term must be 2n — 2 = 2(D + 1 — s) and maximum possible 
rank is iV + 2, therefore we get, 


IV + 2 > i > 2(1V — s + 1). (4.14) 

For this equation to have a solution N < 2s. So we need at max 2s derivative order parity- 
even terms, to construct parity-odd terms till (n — 2 + s) derivative order. For example at 
leading order, s = 1, only pure terms with at max 2 derivatives are required. The parity even 
terms required till s = 3 are enlisted in tables (3) and (4). Further, if we were only interested 
in finding terms that survive at equilibrium, we can use the Killing condition and drop all 
terms with # derivatives. 

Some of the combinations constructed by this procedure using table (3) are: 

1. (2D indices): D( 2,2,1) 

2. (a) (2D - 1 indices): (D - 1)(2, 2,1) ® (1,1,1) 

3. (a) (2D - 1 indices): (D - 2)(2, 2,1) © (§, 3, 2) 

(b) (2D — 2 indices): (D - 2)(2, 2,1) 0 2(1,1,1) 

(c) (2D — 2 indices): (D - 2)(2, 2,1) © (1, 2, 2) 

and so on... The counting can be extended arbitrarily to the derivative order we need. In 
next section we will construct terms till subleading order, and later in appendix (B) we 
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Table 4: Parity-even Data-types - Vanishing at Equilibrium 


Data Type 

Decomposition 

Local Rest Frame 

(2,2,1) 

<?* := 1 SM 

(1,1,1) 

Vj := S i - diD) 

(1,3,2) 

Z ijk := n ijk o - 2 iry k x? + ^ (fi j V k Vo + \X[ k VH 0 - 1 x( k v%) 

^(2,2,1) 

d l a jk 

(1,2,2) 

ST := n' 0 3 0 + IdVXiJ - 2 ^cWtf - ^X{ a x( a 

3q(2,2,1) 

doX{ k , doa'o 

^(1,1,1) 


(1,1,2) 

3 0 (1,1,1) 

9 0 9 0 ^a, <9 0 Vij 

(|,4,3) 

#(§,3,2) 

Qi^jkl 

3q(2,4,2) 

d 0 JC 3klm 

##(2,2,1) 

d l d 3 a kl 

(1,3,3) 

#>(§,3,2) 

d 0 Z i3k 

^(1, 2,2) 

d ir E? k 

#>#(2,2,1) 

d 0 d 3 X ki , d {) d l o 3k 

##(1,1,1) 

d l d 3 Vj 

(1,2,3) 

#>(1,2,2) 

doZV 

#>#>(2,2,1) 

dodoX kl , do do a 13 

#>#(1,1,1) 

dod 3 d k "d\, dod r Vi 

(|,5,4) 

##(§,3,2) 

Qi Qj 77 klm 

#>#(2, 4,2) 

d 0 d j ic klmn 

d i d 3 d k (2, 2,1) 

d l d 3 d k a lrn 

(1,4,4) 

#>#(§,3,2) 

d 0 d l Z 3kl 

d l d ] (1. 2. 2) 

Q i Qj~ j kl 

#>#>(2,4,2) 

dod 0 JC klmn 

do##(2,2,l) 

d 0 d 3 d k X^\ d 0 d 3 d k (j lrn 

### s (i,i,i) 

&d 3 d k Vi 

(1,6,5) 

###'(§,3,2) 

q i Qj Qk^lmn 

# ) ## ! (2,4,2) 

d 0 d 3 d k IC lmno 

d t d 3 d k d l (2, 2,1) 

d l d 3 d k d l <7 mn 


will extend it to subsubleading order. We will suppress the usage of data-type (2,4, 2) for 
brevity; combinations involving it can always be reached by exchanging (2,4,2) with two 
(2,2, l)’s. 




4.3 Examples of Parity-odd Counting 


4.3.1 Leading Order (D=n-1) (s=l) 


For s = 1, the required indices are merely 2D = 2n — 2 (V e ), which amounts to the only 
combination: 

£>(2,2,1), (4.15) 

along with the terms involving (2,4,2). However in V e all the free indices are contracted 
with Levi-Civita, which will kill any term involving (2, 4, 2) due to Bianchi Identity. The only 
remaining combination is - (n) vectors 

/m— 1\® n 

\n—m -i ’ 

m —1 

where we define, 


m 

n — a — m 


m 


1 


m+a 


n y^xVx I I y^yVy 
x=a -\-1 


y=m-\-a -\-1 


l 


m+a 


n _ a _ m/ 2n _ o eiwj2...^„ 

ll 232 •••I'aja x=a-\-l y=m-\-a +1 


n n 


(4.16) 


4.3.2 Subleading Order (D=n, s=2) Surviving at Equilibrium 

At subleading order, index families required are: 2D = 2n (\+), 2D — 1 = 2n — 1 (T e ) and 
2D — 2 = 2n — 2 (V e ). We only compute terms surviving at equilibrium because that is what 
we need for the current work. 

2D Family: 2D family was already discussed in § 4.3.1, but this time since two indices 
are free from e, one (2, 4, 2) can appear with two antisymmetric indices of R l i kl contracted. 
However we are supposed to take a contraction on remaining indices, which again due to 
antisymmetry vanish. Only remaining data are - (n — 1) vectors: 

i • , 71—1 

\n—l—m/ijk 1 2a m=1 

2D-1 Family: Combinations in (2 D — 1) family which survive at equilibrium are: 

1. (D — 1)(2, 2,1) © (1,1,1) 

2. (D — 2)(2, 2,1) © (|, 3, 2) 

3. (D — 3)(2, 2,1) © (|, 5, 3) 

along with the combinations with (2,4, 2). In T c only one index stays free from e, hence again 
(2,4,2) and (|,5,3) cannot appear. The remaining two combinations will yield: 
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1. (n — 1)(2, 2,1) © (1,1,1): 2 possibilities - (6 n — 4) traceless symmetric tensors and (2n) 
scalars 




-1 \(b' fc ; 


0 


n—1 


m=1 


Scalars: 


c-tyw a 


n 

m= 1 


2. (n — 2)(2, 2,1) © (|, 3, 2): 1 possibility - (2n — 2) traceless symmetric tensors 




Ajk 


n—1 


m= 1 


2D-2 Family: Combinations in (2 D — 2) family which survive at equilibrium are: 

1. (D — 2)(2, 2,1) © 2(1,1,1) 

2. (D — 2)(2, 2,1) © (1, 2, 2) 

3. (£>- 3)(2, 2,1) ©(§,3,2) ©(1,1,1) 

4. (-D — 3)(2,2,1) © (|, 4, 3) 

5. (D — 4)(2, 2,1) © 2(|, 3, 2) 

6. (D — 4)(2, 2,1) © (|, 5, 3) © (1,1,1) 

7. (D — 4)(2,2,1) © (|, 6,4) 

8. (D — 5)(2, 2,1) © (|, 3, 2) © (|, 5, 3) 

9. (D — 6)(2, 2,1) © 2(|, 5, 3) 

Along with these, we have the combinations with (2,4, 2). However, V e has no index free from 
e, and hence Bianchi Identity will not allow (2,4,2), (|,5,3) and (|,6,4). Further, (1,2,2), 
(|,3,2) and (|,4,3) will vanish as they cannot be made completely antisymmetric. Finally 
only one combination will remain, yielding: 

1. (n — 2)(2, 2,1) © 2(1,1,1): 1 possibility - (n — 1) vectors 

At equilibrium we have (2n) scalars, (2n —2) vectors and (8n —6) traceless symmetric tensors. 
We have tabulated these data and their equilibrium values in table (5). 


4.4 The Basis of Independent Data 

As we discussed in § 4.2, the data we have enlisted in the preceding sections is a ‘complete 
set’ but not independent. There might exist numerous relations among them through anti- 
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Table 5: Independent Leading and Subleading Order Parity-odd Data at Equilibrium 


Name 

Term 

Equilibrium 

n 

m.=1 

\n—m) Uu 

lm— 1\ 1 
\n—m/ 

Sa?71 

n 

m,= 1 

CVp&A 


vL 

n —1 

m= 1 

\n— 1— ml ^vpcr 1 

/ m-1 \ fiafj 

\n— 1— m/kijJl J 2a 

vL 

n—1 

m= 1 


{XXXXCXdho 

r 

1A m 

n 

m=l 

l { rttP u)a V a # A 

a» 

2 Am 

n—1 

m= 1 


h jt 

3Arm 

n—1 

m= 1 


L-yy^jiAjy 


symmetrizations of 2 n or more indices. If we look back at § 2, the need of all independent 
data arose to write down the most generic form of the constitutive relations. We write the 
energy-momentum tensor and charged current as a combination of all independent tensors 
and vectors respectively up to some undetermined coefficients which are called transport co¬ 
efficients. We then determine the same quantities from equilibrium partition function and 
compare with the fluid results. It turns out that the transport coefficients which destroys 
the positivity of entropy current divergence are set to zero by this procedure. We call these 
transport coefficients unphysical. Put differently, the partition function generates only the 
physical transport coefficients in the constitutive relations (eqn. (3.14)) at equilibrium. 

Now if we relax the condition ‘independence’ while writing fluid constitutive relations, i.e., 
add more terms to these relations which could have been determined in terms of others; they 
can be regarded as redundant transport coefficients in our system. Since the charge current 
and the energy-momentum tensor we derive from the partition function remain unchanged, 
we get relations between the transport coefficients (including the redundant coefficients) and 
the coefficients appearing in partition function. However, we still have our answers - the 
independent transport coefficients and distinct constitutive relations. 

Let us explain with an example. Suppose at some particular derivative order, we have total 
I number of vectors V} 1 . We can write charge current at this order as, = J2i=i a i^t > 
where Oj’s are transport coefficients. On the other hand, suppose our partition function 
has X number of independent coefficients Cj’s, and it generates a charge current = 
LLi *{C s )V?. Ci(Cj) are some functions of Cj’s. By comparison we will get a* = Cj(Cj). 
These are I relations with X free parameters, and thus imposes I — X constraints on a*. 

Now let’s add to our set K more vectors V£, a = I + 1,... ,1 + K which could in principle 
be determined as: Va = J2i=i C a iVh . Then we would have guessed our ansatz to be = 
1 biV^, and by varying partition function we will get = Yln=i diiCj'jVj 1 . di(Cj) are 
some functions of Cj’s determined by relation a = (di — d a C Q ij, as our partition 


- 24 - 


























function is still the same. By comparison we will get = di{Cj). These are K + I relations 
with X free parameters, and thus imposes K +1 — X constraints on 6*. We hence get exactly 
K extra constraints, to kill the K extra degrees of freedom we added in the system. But 
once we have imposed these constraints, we will only be left with X independent transport 
coefficients. 

However, note that we still need independent set of scalars that enters the equilibrium par¬ 
tition function, for our arguments to make sense. We check it here before we proceed. At 
leading order there are no scalars. At subleading order the scalars do not have enough indices 
for 2n or more antisymmetrizations, as a result all the scalars we get are independent. At 
higher order however, it may not be so easy to find out all the independent set of scalars. 

Lets look at an example of such residual 2n-antisymmetrization conditions. In eqn. (4.13) if 
we chose B to be of the form S g lJ , we will get: 

g(_ ir +l e <K..A.- M . .^ = 0. (4.17) 

a=l 

where ( ) denotes the traceless symmetric part of a matrix. Hence one of these matrices of 
type T c (after making traceless) is not independent for a given A. A similar argument is 
valid on other tensors like Tp using S c But as we are treating all symmetric traceless 
tensors (of type [ ] e ) to be independent, this should reflect in our final constraints, and as 
we will see, it will. It turns out that till subleading order, eqn. (4.17) is the only remaining 
residual constraint, and thus we can construct an independent basis; but this issue might turn 
more subtle at higher derivative orders. To illustrate the procedure we will not start with the 
independent basis even for subleading order, and show that we get consistent results at the 
end. 


5 Fluid Constitutive Relations 

Having all the data we require, we are ready to find the constitutive relations for fluid. We 
start with the results which are already known in literature, i.e. fluid up to leading derivative 
order. We revisit the results in our notation. Later we consider charged fluid at subleading 
order in § 5.3. We also set up the notation and architecture for subsubleading order parity- 
odd fluid in this formalism in appendix (B). However we do not compute the constitutive 
relations explicitly, as we will discuss, the calculation becomes a lot non-trivial. 

5.1 Ideal Fluid 

At zero derivative order only energy-momentum gets a transverse contribution: 

= AP^, (5.1) 
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where A is some arbitrary function of •& and v. Now comparing eqn. (2.8) with eqn. (3.14) 
we can write at ideal order, 

E oV y 0 + A 0 y = yp 0 , 

- E 0 V l 0 J 1 + Viovi - V ° + 'doVoQoVo = 0, 

y/l + VioVl 

£ ” (1 + ^ , + (TT^)= e ”' 

Q 0 v l 0 = o, 

Q\Jl + VioVl = q Q . (5.2) 

The identifications will then give 

v 0 = -e CT , v l Q = 0, A = P, Q = q, E = e. (5.3) 

Note that we have identified A, Q, E exactly, and not just at equilibrium, as we explained in 
§ 2. Therefore, the energy-momentum tensor and charge current for ideal fluid can be written 

clS ? 

T$ = evPu v + PPV", = qvP. (5.4) 

5.2 Leading Order Fluid 

One can divide the constitutive relations in hydrodynamics in two different sectors - parity 
even and parity odd. Ideal fluid belongs to the first sector (in d > 2). The first non trivial 
derivative corrections in parity-even sector appears at the first derivative order e.g. shear 
viscosity term in energy-momentum tensor. Whereas in the parity-odd sector, the leading 
terms appear at (n — 1) derivative order for a fluid in 2 n dimensions. All these terms and 
the corresponding transport coefficients (at leading order) have already been found in [6]. We 
shall discuss their result in our notation. 

5.2.1 Parity-odd 

Since there is no parity odd scalar and transverse symmetric traceless tensor at (n — 1) 
derivative order (see table (5)), only charge current gets parity-odd corrections: 

n 

= Y, (5.5) 

m =1 


The combinatorial factor is introduced for convenience. It also ensures we do not surpass the 
limits of m. The fluid variables receives following corrections, 

■&A = d Ao + A (n-1 >0 A , «* = «*+ A (5.6) 
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Further, there is no parity-odd gauge invariant scalar at equilibrium on M. 2 n-ii implying that 
= 0. Now comparing eqn. (2.8) with eqn. (3.14) we will find the constraints at 
parity-odd leading derivative order: 


UJm — 


$ 2 n 
e + P 


sC m —l + QCm + (jl + 1) ( 1 -■—— 

m m + 1 


c (2n)i 


(5.7) 


And the corrections to fluid variables, 

n 

= A (n- l) v = Q) ~^n-l) v i = £ n-l C m _ iao(m) £ (m) 

m— 1 


where, 




il 2 n 

e + P 



W ~ Cm + 


C ( 2 n) ^m+1 

m(m + 1) 


(5.8) 


(5.9) 


Here we present these relations for completion as well as to set up our notations and conven¬ 
tions. We would also like to make some interesting observations about these functions. One 
can verify that 

sa m + qoim+l = V'oJrn ~ ^m +1 V 771 6 {1,71 - 1} , (5.10) 

P (1 ’°W = s (p (1,0) a m ) ( ’ ) + q (p (1 ’ 0) a m +i) ( V m G {1, n — 1} . (5.11) 

Here pressure P($, z/) is function of temperature z? and redehned chemical potential za For 
any function Q($, z/) we define Q( m ’ n ) = -Q. These will come handy in subleading order 

calculation. 


5.2.2 Parity-even 

The most generic current corrections at parity-even leading derivative order are (see ta¬ 
ble (1)): 

3 

T (1) = E A a^a , n (T) = -2^ - C P^e, (5.12) 

A=1 

while at equilibrium the only surviving contributions are: 

+(i) = E A °a+ a , = 0. (5.13) 

A=1 

There are no gauge-invariant parity-even scalars at equilibrium that appear at this order. 
Therefore, A= 0. Now comparing eqn. (2.8) with eqn. (3.14) we will find at 
parity-even leading derivative order that all corrections vanish 

7rjj 1} = ^(i) = A (1) i? = A (1) z/ = A (1 V = 0. (5.14) 
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We hence get the constraints: 


Ai — A 2 — 0. 


(5.15) 


So finally the form of currents is 


r I'j } = A3V3 1 , ngj = -2 rja^ - (p^e. 


(5.16) 


We also get to know that no fluid quantities ( 1 ?, v, v l ) get order one parity-even correc¬ 
tion. 

5.3 Subleading Order Fluid 

In this section, we shall describe the constraints on charged fluid in arbitrary even dimensions 
at subleading derivative order (i.e. n order), in presence of U( 1) anomaly. Where as, the 
subleading correction to parity-even sector comes at second order in derivative expansion. 
Some aspects of four dimensional fluids at sub-leading order have already been performed in 


[7, 16]. 


5.3.1 Parity-odd 


Sub-leading order parity-odd fluid dynamics in four spacetime dimensions has already been 


discussed in [7]. Here, we generalize the results in arbitrary even dimensions and find the 
constraints on the transport coefficients. We see that, much like in [7], the higher dimensional 
transport coefficients depend on first order transport coefficients 77 , £. 

From counting we can see that the n order parity-odd corrections (at eqb.) are given by (see 
table (5)) 



(5.17) 


m= 1 
ri—1 


n —1 



m =1 



o2,AFm 



n 



(5.18) 


Sum over the relevant ‘A, T’ indices is understood. We explicitly write the m index contraction 
to emphasize that the sum runs over different values for different terms. We do not state non¬ 
equilibrium contributions as they won’t be required in this computation. 

From eqn. (2.8) and eqn. (3.14) we get, 



(5.19) 
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AW J = -e a A^q = —€ 


<5w; 


egfe 


(C) 




(5.20) 


Now, 


A( " >£= (^) o Al “ , ’ ,+ Oo A< " > ' / ’ Al " )<,= (i) o AI “ ),,+ (^) o A< " V <5 ' 21) 


Therefore from eqn. (5.19) and eqn. (5.20) we can write, 

A (n R = ^ 0 E oAr .. (C) , 
Rr 


where, 


/.□Ml iwi ^ 

pi _ [ 17 de I qr’ i? <9g le 

^ Ar “ ,4 9^ I 1 I 

\ 17 de I q ’ i? 3g I e , 


= R^} • 


(5.22) 


(5.23) 


Similarly comparing the T l equations in (2.8) and (3.14) we get, 

° + (C) (c) 




SA; 


— e 


<9a, ; 


which can be written as, 


where, 


, . . u A S&Wwfg 

AW«' = (-) -^ (c) 


p v 

^ o 


(bo) ^a A - 


(5.24) 


(5.25) 


A a = drPE rA = (tff | 5> | e ) , hA = R^} , a\ = {da\ A*} . (5.26) 


One can check that E A r is symmetric matrix and d\ = . We would like to emphasize that 

these are purely notations, to make the calculations tractable and easy to digest. There is a 
summation on repeated A, T indices. Now comparing T lJ and J l in eqn. (2.8) with eqn. (3.14) 
at parity-odd subleading derivative order, we have corrections to constitutive relations 


#KU~ 2 - 


5 A Ww^ 


P?' 0) ii(n) = VoSoA 


&9ij 

5A^Wft 


SAWwffi i .. .. 

- g lJ A oA f Q (C) - -A^tt”, 


6ci\i 


S'doA 


(5.27) 


where, 


o _ dP 
0 

<AA 




(5.28) 
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A swaps the value of A : 1 2. The generating functional A^Wf^ contain all scalars S oAm . 

But one can check that S 0 i m can be connected to S 0 2 m by a total derivative. So we take the 
partition function 

/ n 

d 2n ~ l x^Y n_1 C m -i5 m S o2 , m . (5.29) 

m=1 


We compute the variation of generating functional with respect to different fields and find 
that 


8 A 

$9ij 

8 A 

S r d 0 A 
8 A 

8a m 


= 0, 

n 

= -(-) A E n - 1 C m -l4 1 ’ 0) SaA, m> 

m= 1 
n —1 

= (»- 1 ) E AATavL 

m= 1 


(5.30) 


Using the form of lower order currents corrections from eqn. (5.16) we can write, 


= -2r] 0 A^-^a ij - C 0 g ij A^-^Q 

= —2rj 0 i} 0 ^Cm-tdA (^) Tjj Am - Vo n ~ 2 C m _!(n - l)a 0(ro+2 _ A) T^ iAm 

- 9 ij (o#o "-'Crr^ldA S oAm (5.31) 

A (n_1) <j(\) = X o3 A^£ i 

^o3 ^m— l(p* 1) [^o(m+l) "h ^olm’ (5.32) 


One can now use the results, obtained in eqn. (5.30) and eqn. (5.32) in eqn. (5.27) and 
comparing these expressions with eqn. (5.18) to get the constraints, 


Tl,A m — ‘2T) r &d\ ^ J ) ^2,Am — TJ^Tl l)^(m+2—A)? ^3,Arm — 0, 

■ Aa ~ l/ fl ’ lm ) = at (^ 2m ~ i/ fl ’ 2m ) = p(1,0) * m ’ 

hm = —A 3 (n - 1 ) (a m+ i + z/a m ), z> 2 m = -(n- 1 ) (g,\ m +i + s% m ) • 


(5.33) 

(5.34) 

(5.35) 


Hence everything is determined in terms of a known function a m and a new coefficient \ m . 
Note that if we had used the 2n-assymetrization condition eqn. (4.17) to get rid of one traceless 
symmetric tensor to start with; a consistent choice would have been to remove T 3 Alm entirely 
and T 3 A2m for m = 1 (see table (5)). The coefficients of these terms are set to zero already 
by our constraints, which means the other leftover constraints are independent. 
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Finally we get the corrections to fluid variables using eqn. (5.25) as 

A (n) 4 = n_1 C m _ 1 (—) r E oAr x' o(m) S 0 p 

n —1 

A (n V = (n - 1) n ~ 2 C m -i (xo(m+ 1 ) - z'oXom) vl 2) m- (5-36) 

m =1 


5.3.2 Parity-even 

Next, we present the results for sub-leading order (two-derivative) parity even sector for the 
fluid. From counting we can verify that at the second order, parity-even corrections (at eqb.) 
are given by (see table (2)): 

rS ( 2 ) = = + < 5 - 37 > 

# ° # # 


# refers to sum over all relevant indices. Now comparing eqn. (2.8) with eqn. (3.14) at parity- 
even subleading derivative order, and performing a similar manipulation as last section, we 
have corrections to constitutive relations: 


1 ij n 

°(2) - 2 “ 


6AWW$ 


8f_h 


-g ij A oA - 


6A&W$ 




<5$oA 


p( 


+ 


g ij A^v k |A o2 ^ (1) + (i-A^Pj 1 - 0 ) + ^-A^Pj 0 - 1 ) - ipi 1 ’ 0 )) A(V j 


^ (1 ’ 0) 4( 2) = ^oSoA 


(C) pa.olA^)^) 


5a\i 


(5.38) 


while the fluid variables get the corrections: 

(5A( 2) 1F^ 


A (2) d A = ^ 0 E oAr - 


(C) 


(5-dor 


a„a - A^Ai'V 


doEoA 2 A (1) u4 (1 ) 


r O 


(bo) 5 a m 


(5.39) 


Notice that the boxed terms only contribute for four dimensional fluids (n = 2). Out of 
the scalars enlisted in table (2), S 0 i A can be related to others by a total derivative. Hence 
is given by: 

= J {dx‘}v5{5«S„ 4 + ‘5/ArS o3 ( A r) + Stf A rS 0 2(Ar)} ■ (5.40) 
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Now we can find the variations of A^W^y 


8 A 


( c) 


tig. 




= -9a5rT^ 1a - (c?Ai9r Sr - S$at) T o 2 ,Ar + 2 ^/ArT ^ 3 Ar + SrT 


o4 


+ g 
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1 — ^ ^ OaSrSc^a + (^d a dr Sr — ^——j-<9A<9r Sr — -S$Ar + ^ __ y <SflAr ) S o2! Ar 


1 f 4 

_ 2 l 1 ~~ d- 1 


5/ArS 0 3,Ar — y f 1 — d - 1 ^ 


(5.41) 


M 


oE 


SA^W^ / ! \ 1 1 

— <S0eaS o i,a + ( ^(r^i?A)s — ydsStfAr J ® o2 . Ar — 2^ s ^Ar^o3,Ar — - Q^SrSoa, 

(5.42) 

(5.43) 




SciAi 


— 2Sf t \Y V* ir — 2cK]‘S/ArV) ) 2 jrE . 


Using the form of lower order corrections from eqn. (5.5) for n = 2 we can write, 
A^^i) = WoE^oS 

= (—) S W 0 E {MoE a A^ol,A — (/^oE^a) ^o2,Ar} • 


(5.44) 


We can now put the variations of generating functional along with lower order corrections 
worked out above in eqn. (5.38). Using eqns. (5.10) and (5.11) and eliminating partition 
function coefficients S's we will find following 7 constraints, 


7i,A + <9 a t 4 = 0 , 


d -2 


°ol,A = J-7<9 aT4 ~ A^5s3aT4 - AsT2 EA) 

d— 1 


d 3 

2o2,Ar = 3 a^t74 - Ae<9e<9a<9tT 4 - ^ _ ^ r 2,Ar — 2As<9(A r 2,r)s + AsdsT^Ar, 

d — 5 

4(J3 Ar = — - 7 — 7 r 3 ,Ar + A S (9 s t 3 aia 

d — 1 

d- 3 4 „ 

2 cr 4 = — -— 7 X 4 + A s dsr 4 , 


d- 1 

d 

^’’A ~ p(l, 0 ) 

" 2 ’ rA - - p( li0 ) 


SrT3,Ar, 


Ss(9A73,rE- 


(5.45) 

(5.46) 

(5.47) 

(5.48) 

(5.49) 

(5.50) 

(5.51) 
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Coincidently none of the constraints depend on n = 2 special contributions. On the other 
hand fluid variables corrections are given by eqn. (5.39): 


= ??oEoOS 


(t 0 2,ea + <9a<9st 0 4) S 0 i,a + ^<9(r r o2,A)s - -dy,T 0 2 : \r + -<?£<9 A <9rTo 4 ^ S o2) Ar 


-T5sT o3i ArS o3 (Ar) - -<9 s r 04 S o4 


+ 


^oEoosSs (E o (1 ’ 0) a oA a of ) - A on a oA a oT + i? 0 E oQ2 a oA ( \q 0 ot 0 f - w of 


5 o3(Ar) 


(5.52) 


A< 2 V = ( _ } A JfoA 


p y 

r o 


( 1 . 0 ) 


7b3,Ar 


p(l.O) 

-1 n 


a oK a oT 


V 


oi,r 


— d s ( r o3)A r - 


p(fiO) 


a oA a or 


V, 


o 2 ,rs 


(5.53) 


This completes our calculation of subsubleading derivative order fluid. 


6 Conclusions 

In this paper we computed the energy momentum tensor and charge current for a fluid system 
in 2 n dimensions with U (1) anomaly up to subleading order in derivative expansion (for both 
parity odd and parity even sectors) from the equilibrium partition function of the fluid. We 
described a novel counting prescription to construct the fluid data. However, an important 
issue we encountered here is that it is non-trivial to find independent vectors and tensors at 
arbitrary derivative order. But we were still able to find the independent transport coefficients 
and distinct constitutive relations. We showed that the knowledge of independent scalars at 
the required derivative order is sufficient for this purpose. This is a powerful observation and 
it enables us to carry on the computation at (n + 1) derivative order, where, we could find 
the independent scalars. We observe that the parity odd transport coefficients which appear 
at n derivative order in constitutive relations are constrained and some of them depend on 
the first order transport coefficients like rj, ( etc. It would be interesting to find the similar 
dependence in a holographic set up [17]. We plan to explore the holographic computation in 
future. 

It is also interesting to find the fluid constitutive relations in presence of both U( 1) and 
gravitational anomaly in arbitrary 2 n dimensions. But, since the gravitational anomaly ap¬ 
pears at two higher derivative level compared to the (7(1) anomaly, it requires to carry on 
our analysis to one higher derivative (sub-sub-leading) order, i.e. to (n + 1) derivative order. 
Fortunately, as mentioned earlier, even at this order, we could determine the independent 
scalars and hence, in principle, the computation is possible. We have carried a large part of 
it in appendix (B). 
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A Kaluza-Klein Decomposition 


If a (d+l)-dim spacetime Aiu+ 1 ) has a preferred time-like direction uU, it can be decomposed 
into S 1 x Awhere S' 1 is the euclidean time circle. A fc-rank tensor decompose in 2 k parts 
in this scheme: 


1. S does not decompose. 

2. decompose in oo^V^ and P^ U V U . 

3. T^ decompose in A‘P™T ]1V , PT^ v T ]1v and P^P^T^, 

and so on. Where Po V = G liV — c ■jp is the projection operator. If we are studying theory 
at equilibrium, we already have a preferred direction along the Killing vector of the theory 
uj ,j = do. In this case we know that a (d + l)-vector will yield a scalar: 


=> Vo := V, 


(A.l) 


and a (d)-vector: 


P^V U => V := V\ 


(A.2) 


Hence we see that a U( 1) gauge field A M will be decomposed in |Ao(x), A*(x)}. Similarly a 
tensor T^ v decomposes in Too , T* 0 , T 0 l , T l P It is the similar way the metric G^ v on 
decomposes, hence we define: 


G 00 = -e 2CT , G l 0 = 0, G ij = <f\ 


(A.3) 


where we define g lJ as metric on M^y Now using the diffeomorphic invariance one can work 
out the full form of G^ 


ds 2 = G /lu dx k 'dx l/ = —e 2a G) (dt _|_ a,i(x)dx l )~ + gij(x) dxMx J , 


= 


2(7 


—e 


—e 2<7 ai 


-e 2a aj 

(. 9ij - aiCije 2a ) 


G IW = 


(—e 2<t + a 2 ) —a J 


—a 




(A.4) 

(A.5) 
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where time redefinition invariance requires that a 1 is an independent gauge field, named 
as Kaluza-Klein gauge field. Using the euclidean time period ft we can define the local 
equilibrium temperature of the theory as: d 0 = l/ft 0 = e~ a /ft. Our higher dimensional 
metric is hence disintegrated in a scalar (Temperature), a gauge field and a lower dimensional 
metric. 

We can now use the metric G^ u to raise/lower the components of vectors: 

H = gij V j + a,;Vo, V° = -e" 2 ^ - aj V j . (A.6) 

which are not Kaluza-Klein gauge invariant. From here we read out the (d)-covectors: 

U = (V, - a*V 0 ). (A.7) 


Determinant of metric in two spaces can be related as: 

G = — det G^ = e 2<T det g l3 = e 2u g. 

We have the Levi-Civita symbol in lower spatial dimensions: 

e ijk - = e a e 0ijk - = -e~ a e 0 ijk -, 
where e 0123 ’" = \j\[G and e 123 '" = 1 /^/g. 

It is useful to see how higher dimensional contractions behave in lower dimensions: 

= —e~ 2a AB + A l B t 


(A-8) 


(A-9) 




A 




— e -V e lJl32-3n-2 


l Y A jl-ia-lQja-. 


•Jn — 2 * 


(A.10) 

(A.ll) 


A.l Derivatives of Metric 


Once the metric is known we can reduce the derivatives of metric, i.e. the Christoffel Symbol 
and the Riemann Tensor. The Christoffel Symbol is defined by: 



G Xa {d^G au + d„G afi - d a G^) 


(A.12) 


Pretending it to be a tensor at the moment, if we define its indices to be raised and lowered 
with the metric G^ v . We can reduce it for Kaluza-Klein form of the metric: 


Vooo — 0, 


V'oo = -e 




2a u u ° 


do 


P i _ ■p i _ ^2cr u u ° 

1 0 0 ~ 1 00 — e ^ 


r y 0 = 0 j = -e 2<T /* J , 


V J ' = ~ 2 e2 VV 6 ( d a a b + d b a a ) 


f kij =g il g jm T k ltn , (A.13) 
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where T k - is Christoffel Symbol on which is raised and lowered by g l3 . Also we define 

KK field tensor: 

/v = VV-W. (A.14) 

r 0 ' J is not KK gauge invariant, even though it has time index down and spatial index up, 
which is the manifestation of T not being a tensor. 

Lets define the higher dimensional covariant derivative as V and lower dimensional as V, 
whereas the usual derivative is given by d. We can check that: 

VA = VY J + 

VWo = VV + p o VV% + 

2 

V°v* = + \fFpVj, 

V 0 V 0 = e 2CT W^, (A. 15) 

Vo 

similarly, 

yiyjk = yiyjk + I/b'v o fc + I f ik V j 0 , 

V’V J 0 = + \f ij V oo+ = + VavA 

VWoo = VWoo + 2V 00 ^ + \/ a f \(V + V A 0 ) , 
v °yd = + -e 2 "AV ife + ^ W’ 

V 0 V l 0 = Voo^ + ^e 2CT AV fc 0 + e 2 ^^, 

VoVoo = e 2CT (V* 0 + V 0 l ) (A. 16) 

VO 

Finally the Riemann Curvature Tensor is defined using an arbitrary vector X 11 as: 

K^paX” =^(vA-VA)l„ (A.17) 

using which we can define: 

A- = A-A n = n a a . (a. is) 
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Now a straight away computation will yield: 

n = R- 44 VitfoVMo + 2-j-V i 'Vi'& 0 + 

'0 2 o i3o 4 

= e- 2CT ^ 00 = 24 v^ 0 V ^ 0 - 4 ViV^o + 4 2,T f%-, 

"do 4 

= e-^'o = e CT i (v fc / fci + , 

^ - 24v^oV^o + 4Wtf 0 + 

vi Vo 2 

« a ^^i ^ = e-^o'o = 2 iw, - + \e 2 °f a r, 

n^ ka u a = e-°K ijk 0 = ^V k fl j - 4 (/r'V^o + ^/{Vi ? 0 - i/f fc V^ 0 ) • (A.19) 

Here is defined to be lower dimensional Riemann tensor, and HA = B ll j ,, i? = R i i . 

A.2 Derivatives of Gauge Field 

Now let us have a look at derivatives of gauge field A ,J; . Being a vector it decomposes as: 

A = Ao = -e 2a (A 0 + aj A>) , A i = A \ Aj = (Aj - ajA ). (A.20) 

The gauge transformation A M —> A^ + d^A translates to: 

A —> A, Aj —> Aj + djA. (A.21) 

Hence A 1 is a gauge field on Jv[ ( c n , while A is a scalar. Using /I (euclidean temperature) 
we define the local equilibrium potential v Q = /3A. Higher dimensional field tensor however 
decomposes as: 

jrvv = ^ jrij = pij + e °$ oUo fij^ jri Q = e CT ^ 0 VV 0 , (A.22) 

where, 

V',1' v-'.r. (A.23) 

Now we define the four vector electric field: 

= T^ Uv => £ 0 = -e CT i? 0 u ? ;VV 0 , £* = -e- 2ff ue ff i? 0 VV 0 + Vj (i™ + e a $ 0 v 0 f ij ) . 

(A.24) 
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B Subsubleading Order Fluid 


In this appendix we extend the counting discussed in § 4 to subsubleading order fluid. We 
form a complete set of data at this order and classify the respective scalars, vectors and 
symmetric tensors. Later using the independent scalars at this order we construct an equi¬ 
librium partition function and compute its variation. We were however unable to process the 
constraints explicitly, as the calculations are analytically intractable. 


B.l Counting at Equilibrium 

At subsubleading order (D = n + 1, s = 3), index families required are: 2D = 2n + 2 (V^ ), 

2D - 1 = 2n + 1 (Tf), 2D - 2 = 2n (Vf), 2D - 3 = 2n - 1 (T e ) and 2D - 4 = 2n - 2 
(V e ). We only compute terms surviving at equilibrium, as non-equilibrium pieces are not 
required till subsubsubleading order parity-even or subsubsubsubleading order parity-odd 
calculation. 


2D Family: 2D family was already discussed in § 4.3.1, but this time since four indices are 
free from e, two (2,4,2) can appear with two antisymmetric indices of R^ kl contracted. We 
will find 3 combinations - (19n — 20) vectors of type Vy . 


1. 2(2,4, 2) 0 (n — 3)(2, 2,1): 1 possibility - (n — 2) vectors 


/ m— 1 \ijklm abv~ 

\n— 771 — 2 / 'S'fc ^Imab 


n —2 


m— 1 


2. (2,4, 2) © (n — 1) (2, 2,1) - 4 possibilities - ( 8 n — 10) vectors 


/m—l\ l y " / m— 1 -y Y'a / rn—1 \ l 3 K y ab y 

\n—m) ^ , ’ \n— 771 - 1 / ^Aja^ \ /V A\ ab 

m= 1 tti=1 J 

/ 777—1 \ijklm ab y y 
\n-m- 2/ ^jk XhalAVbm 

3. (n + 1)(2, 2,1) - 3 possibilities - (lOn — 8 ) vectors 


77—1 


-1 ©L 

;-l/ 
n —2 

tti=1 


n—1 

771=1 


/77i—1\ 2 W/c -y / 771—1 / 771—1 \ijklm - , -pfr 

\ 71 — 771 / ’ \ti— 771 —1/ J \n— 771 —2/ 


771=1 


71—1 


m—1 \ijklm . 


771=1 


71—2 


771=1 


2D-1 Family: 2D — 1 family was already discussed in § 4.3.2, but this time three indices 
are free from e. So only one among ( 2 ,4, 2 ) and (|, 5,3) can appear, and not more that once. 
We will find 5 combinations of type : 

1. (2,4,2)®(n — 2)(2,2,1)©(1,1,1): 3 possibilities - ( 8 n —12) symmetric traceless tensors 

771—1 
71—771—2 


a) 


71—1 


771=1 


-1 \<bfc^a>6^ A 


/ m—1 \ 

\7l—771— 1/ 




71 — 2 
771=1 


71—1 

5 

771=1 
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2. (ra)(2, 2,1) © (1,1,1): 5 possibilities 


(a) Contraction between (2, 2,1) and (1,1,1) - (12n — 8) symmetric traceless tensors 
and (4n) scalars 

{nZn-f* 


m=1 


m= 1 


Scalars: We can take trace and get 4 n scalars: 

a,** 


771=1 


(b) Contraction between (2, 2,1) and (2, 2,1) - (14n — 18) traceless symmetric tensors 
and (2 n — 2) scalars 


771=1 

Scalars: Taking trace we get 2n — 2 scalars: 

LZ-\) m x 1 ‘x !at d l v A 


71—1 

5 

771=1 


771=1 


71—1 

771=1 


3. (2,4, 2) © (n — 3)(2, 2,1) © (§, 3, 2): 1 possibility - (2n — 4) traceless symmetric tensors 




71—2 

771=1 


4. (?z — 1)(2, 2, !)©(§, 3, 2): 7 possibilities 


(a) Contraction within (|, 3,2) - (6n —4) traceless symmetric tensors and (2n) scalars 


[ n=1 , c-yrxifl 4 *™ 

Scalars: Taking trace we get 2 n scalars: 


i-i \ (ijk y a) ob 


n~l 

771=1 


C--'S9 k X,, u 


771=1 


(b) Contraction between (2, 2,1) and (§, 3, 2) - (20n — 28) traceless symmetric tensors 
and (4ro — 4) scalars. 


/ m— 1 \(ijk \>a)b o -y 
\n—m—l) X A Ob *Tjfc 

Scalars: Taking trace we get 4 n — 4 scalars: 

i-i 


771=1 
71—1 


771=1 


/ 771—1 \{ijklm y,a) 06 v 

’ \ ti — 771 — 2 / ^Aj^kbO 


771=1 

Him 


71—2 

771=1 


771-1 \ l J K y psb y 

i-m-1/ x rjk 


71—1 


771=1 


(c) Contraction between (2, 2,1) and (2, 2,1) - (2n — 4) traceless symmetric tensors 
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n —2 




m= 1 


5. (n — 2)(2, 2,1) 0 (|, 5,3): 3 possibilities - (4n — 6 ) traceless symmetric tensors 


/ m—1 \(b' fc o n 1 / m—1 jr a > ” 1 / m-1 a )b v 

(n-m-l) ^ ifc ©\n-m-l) d ^k U.-2 

J m =1 - 1 


71—1 

771=1 


77—2 

777=1 


2D-2 Family: 2D — 2 family was already discussed in § 4.3.2. Here again, one among 
(2,4,2) and (|,5,3) can appear, and not more that once. We will find 7 combinations - 
(39n — 46) vectors of type Vp: 

1. (2,4, 2) 0 (Z7 — 4)(2, 2,1) 0 2(1,1,1): No combinations possible 

2. (n — 1)(2, 2,1) 0 2(1,1,1): 3 possibilities - (12n — 10) vectors 

„a ak.a U / m-1 \V k 


., C‘i>” ‘XH'P#rMz 

777=1 

m-1 




77—2 

777=1 


77—1 

777=1 


3. (n — 1)(2, 2,1) 0 (1, 2, 2): 2 possibilities - ( 6 n — 4) vectors 


C- 0 'a^A ” ,, (0i^X-,.Wr 


777=1 


77—1 


777=1 


4. (n — 2)(2, 2,1) 0 (|, 3, 2) 0 (1,1,1): 3 possibilities - (16n — 24) vectors 


777=1 


77—1 


%-l \ijk. 


77—1 


777=1 


777=1 


5. (n — 2)(2,2,1) 0 (|,4,3): 1 possibility - (2n — 2) vectors 




77—1 

777=1 


6 . (n — 3)(2, 2,1) 0 (|, 5,3) 0 (1,1,1): No possiblities 

7. (n — 3)(2, 2,1) © 2(|, 3, 2): 1 possibility - (3n — 6 ) vectors 


L-^J ,klm 3aK^X nr 


77—2 

777=1 


2D-3 Family: We are interested in combinations in (2D — 3) family which survive at equi¬ 
librium. We generated them through a Mathematica code and found 22 of them. We won’t 
list all of them here, because it won’t be required. Due to properties of T c , most of them 
will not contribute. We will be only left with 3 combinations - (7 n — 9) symmetric traceless 
tensors: 

1. (n — 2)(2, 2,1) 0 3(1,1,1): 1 possibility - (2 n — 2) symmetric traceless tensors 
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Table 6: Subsubleading Order Parity-odd Vectors at Equilibrium 


Name 

Term 

Equilibrium 

vi l Am 

n 

m ,=1 

^Sia 

CV>'v‘V t tfA„ 

^2(Ar)m 

n 

m— 1 

^S 2 (at) 


v 3(Ar)m 

n 

m= 1 

^S 3 (Ar) 

/m-l\ l nab f 

\77—777/ ■'A «/ra6 

vL 

n 

m= 1 

PS 4 


^5Ar m 

n—1 

m= 1 

/ m— 1 \/“'P°' 

\n-m- 1/ ^i/VApVira 


^6Ar Sm 

77 —1 

777=1 

/ m-1 

\n—m—1/ ^^ApVsrSo- 

(„”©r‘ v J'Wr fe >V“^o 

^7ATm 

77 —1 

m= 1 



^8Am 

n—1 

771=1 


/ \^ k f, pa 

\ 71—1—771/ /Aja-ft fc 

v 9(Ar)Sm 

77— 1 

777=1 

/ m— 1 \fj-vpcr \>ct/3 

\n-m- 1/ Uv^-Apa^-TcrP^s 

/ m-1 \V fc r p fab 

JAjaJrkbJj; 

v lOArm 

77—1 

777=1 

L--- I r l ’°^v l ;v a x rpa 

CTTV’VVVTT 

Vi 1Am 

77—1 

777=1 


VTvbTVVTir 

r 12Am 

77—1 

777=1 

/ m-1 \P^po- a /3 

\n—m—l) u vA'pv *A a/3 

/ m-1 \*l fc o ab f 
\ ?i— m-1/ n ik /Aab 

Via m 

77—2 

777=1 

< n ! n “l 2 > MW ^^i K p^V Q i9 1 V^ 2 


Vi 4m 

77—2 

777=1 

/ "i— 1 \pvpaap T^K y 

\ n - m -2/ ^u^-lKp^-2a /l lSa‘^20 

/ 777—1 \ijklm r na -f fb 

\77—777—2/ l a j J 2k J ^blj2m 

v 15ArSm 

77—2 

777, = 1 


( n !"-i 2 > ufcim /Aa 4 v“/™v m i? So 

Vl6(Ar)m 

77—2 

777=1 


<„”V 2 >‘ J “’"W/AjiV“/n m 

Vl7(Ar)m 

77—2 

777=1 

/ m— 1 \Hvpcrap V K( 5 , , V 

\n-m- 2/ u v^po ^Ana^VSp 

/ m—1 \ z 3 klm a j ) 

\ 7 i—m—2/ - fl 7'fc J AalJTbm 

V 18m 

77—2 

777=1 

7 m7-l \pvpcra0 “ s 

\n— rri— 2/ U u K, p £ K 1C a ps K 

/ TTi-l vijfc/m ab „ 

\n—m—2/ ^jk n Irnab 


Mm!' 

m= 1 

2. (n — 2)(2, 2,1) © (1,1,1) © (1, 2, 2): 1 possibility - (4 n — 4) symmetric traceless tensors 

3. (n — 3)(2, 2,1) © (|, 3, 2) © 2(1,1,1): 1 possibility - (2n —4) symmetric traceless tensors 

2D-4 Family: There are 51 combinations in (2 D — 4) family which survive at equilibrium. 
However none of them will contribute due to properties of V e . 

All the subsubleading parity-odd data surviving at equilibrium has been summarized in ta¬ 
bles (6) to (8). 
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Table 7: Subsubleading Order Parity-odd Symmetric Traceless Tensors at Equilibrium 


Name 

Term 

Equilibrium 

11 1A m 

n 

m ,=1 

tm V 1A 

C-'A’v*/? 

T ^ 
R 2 Am 

n 

m ,=1 

l{P\T v ) 

lm V 2Ar 

C-m> ( ‘/fv t 0r 

11 3A m 

n— 1 

m= 1 

vM 

<„-A, > to ‘/,“/ 2 „W'>tfA 

11 4A m 

n— 1 

m= 1 



^ 5Arm 

n— 1 

777=1 



^ 6A m 

n—1 

ra=l 

/ m—1 \ {^vpa T r ix- a) 

\n—m—l) U "VA pJCa 


R 7 ATm 

n—1 

777=1 

/ rn-1 \^ u P a y a )\r „ 

\n—m—1/ ^ Ap ITcr 

/ m— 1 \ (ij f a) yyb f 

\n—m—1/ A 7 V -frfefc 

11 8ArSm 

n—1 

777=1 

\ra—m— 1/ ^'<A 0 V 2rS ( 7 

C-AANlS/mvMs,, 

9ArSm 

n—1 

ra=l 



u 10 Arm 

n—1 

777=1 


<„ 2 -i 1 > tot /A J W»/f l 

^ 11 Arm 

n—1 

m =1 


/ m — 1 \(b fc A rjbfo) 

\n—m—1/ 4Aj6V 

f'"' 
u 12 Arm 

n—1 

ra=l 

/ m-l \{P v P a „, v a ) S Y7 y 
\n-m- 1/ V^Arpo- 

/ m—1 \{ijk ra)by 7 e 
\n-m-l) * A V ^>/rj/c 

13Am 

n—1 

777=1 

/ m-l \(A“'po’ 7 a ir k>Q 

(n-m-U U u V Aa IC pa 


11 14 m 

n—1 

777=1 

/ m-l \{P u P a - „ k)q 

(n-m-l) U u V a IC pa 

/ m-l p “> 6 

\n—m —1/ V ^ rt 7fc 

11 15 m 

n—1 

777=1 



T ^ 

^ 16Am 

n— 2 

777=1 


<„”-i 2 ) l ’ : ’'“ ra v^ 1 „v fc tf 2o v“i/A Im 

R l7Am 

n—2 

m =1 

/ rri-1 \(pvpcraf3 V 7 n) y. I 

\n—m—2/ «i/A lp A25 cr V Tyafl 

<„-A 2 > < ‘ lMm /i,Ai.w“>/Ai m 

1 18Ar)n 

n—2 

m= 1 

/ m-l \{pvpaotfi t7 

\n-m-2/ u vX lp X25<rX Aa V Tp 

/ m— 1 \{ijklm jt f a )\7 q 

\n—m—2/ hjJ2bkJ Ai^ mVTo 

11 19ArSm 

n—2 

m= 1 

/ m-l \ (p^po-a/3 ^ 

\n-m-2/ Uv^ApAraV S^XaP 

j 771—1 \{'ijkl rn nO') n\) Y7 -P 

\n-m- 2/ JKjJTk V bJZlm 

11 20Arm 

n—2 

m =1 

/ m-l \(pvpvap ir y y K ) S 

\n-m- 2/ U u V\ p XY rT SA ' a 0 

/ m— 1 \{ijklm— q f jd a )b 

\n-m-2/ * jtr\ 0 JT kb-K'irri 

Ji 21A m 

n—2 

m =1 

/ m-l \ (pvpaaP ~ k)S 

\n-m- 2/ UvV 8 ?CA P cTA' a p 

/ m-l \{ijklm „ a>6 

\n-m-2/ ^bjAjk^im 

Ji 22Am 

n—2 

m =1 

/ m-l \ (p^pcct/3 - ^ k ><5 

\n-m- 2/ U u A.AS p VaA, a g 

/ m—1 \{ijklm„ — „ a)6 

\77—777— 2/ J^bj k l m 


Table 8: Subsubleading Order Parity-odd Scalars at Equilibrium 


siA-ir, 1 , 

KnVA, 

{™-l^haifZiV k VAo 

§ 2 Arm 

/ 771—1 \P l/ P cr „ l y yjS y 

\n—m—l) u v<X-A fiSV ^Tpa 

b-Aif ‘AfV^ 

§3Am. | m=1 

Cv 1Ap 

C-™>’v‘/a« 

§4Arm| m= i 

CV2AFp 

<” _ ™)'/A«V^ro 


Independent Scalars 

As we discussed in § 4.4, we only need to construct independent scalars which enter in equi¬ 
librium partition function. At subsubleading order one can find antisymmetrizations which 
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will determine § 3 A m and § 4 Arm in terms of §iAm and § 2 Arm respectively: 




m=l 


= 0 , 


y[hjl •y^m— 1 jm — 1 v^mjri 

cLy ... a >2 


■vs^n — ljn — l 

A. 2 


v fe ^ a]in 


= 0 . 


m= 1 


(B.l) 

(B.2) 


Each of Si Am and § 2 Arm, on the other hand is a unique scalar per choice of the parity-even 
tensor used to construct it by contracting with e. Since antisymmetrization conditions cannot 
alter the tensor structure, these scalars are independent. 


B.2 Attempt for Fluid Constraints 


In the equilibrium partition function A^^Wf^s we can include the scalars: SoiAmj §o 2 Arm- 
But it can be checked that antisymmetric part of S o 2 [Ar]m can be related through a total 
derivative to § 0 iAm- So we have: 

A (n+V >W^ = j {x 4 } ^ n - 2 C m _! {Qi,AmS 0 l,Am + Q 2 ,ArmS 0 2 ,(Ar)m} • (B.3) 

Sum over relevant indices is understood. Varying the partition function we will get: 


(C) _ n—2 


5g, 


Cm —1 


U 


Cl,Am (’I’o9,12Am + ^o8,21Am ) + 2Q 2j ArmT 0 J 7 Arm 


n— 3/ 


C m -l(n 2) Ql,AmT 0 J 182Am Q2,Ar(m+2-E)’S’ 0 ‘l9 ,ASrmj ’ 


-azr 1 = °~ 1 Ql,Am§o2,[12]m d - ^AQ2,rEm§o2,(TE)m J ) 

(B.5) 


SA^^Wf# 

(C) _ n—2 


fiCLQi 


Cm—1 


( fdrQl,AmT o6jrnAm ( ) iZ Cl,Am ( V^ nAm + V 4 ol0jAf7m 


1 . 


— 2dzdrQ2,nAmy l o6,XArm + drQ2,QAm (^V 4 o5 rAm + 2V l o7 AFm + VhofAm 

+2Qo t nAm (2 Vq 8 ,A m d~ ^"olgAm d"~ ^ 012 ,Am 

+ (n — 2) 3 C m -i 


-2 ( 5rQ2,nA(m+2-E) - ^9 r Q 2 ,AE(m+2-n) ) V ol5,(EA)r 


d - f C2,QA(m+2—E) ^ C2,AE(m+2—Q) J ^^ol6,(EA)m A^ol 7 ,(EA)m^ 


(B. 6 ) 


On the other hand from counting we can see that the third order parity-odd corrections (at 
eqb.) are given by (see tables ( 6 ) to ( 8 )): 


E n r = E + p r E r»#s^. (b.t) 


r ^'o(?i+l) / j ' Po H z " o#) J " L o(n-)-1) 

# # # 
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# corresponds to all the relevant indices. Similar to subleading order, here also we will have 
special contributions for n = 2, as 3 leading order (n — 1) parity odd corrections can combine 
to give a 3n — 3 order parity-odd corrections, which will be equal to n + 1 only at n = 2 
(Remember we are not considering n = 1 case). Now comparing eqn. (2.8) with eqn. (3.14) at 
parity-odd subsubleading derivative order, we have corrections to constitutive relations: 


/r °(«+ 1 ) 


SA^+^W^ 


&9ij 

+ ^A o2 A("-% i U (2) + 


SA^+^W^ o 

(C) + g l ^A oK d k P 0 A^v k A^v k 



<5$oA 

( 2 ) + 

a< 1 ) « 1 d 


- g ij A o2 A^v k (q 0 A^v k - ^ (n _ 1) ) 


^ ( 1 ’ 0 ) C(n+ 1 ) = ^S 0 A- 


- 2Pi 1 ’°)A (r '- 1 ^ (i - -i- (A^Trgj + A^tt^) , 

1 /q \ / 

<5A( n+1 )VF^ r i 

Pi^AWg-^P^AWfe + P) 


(C) 


6<l\i 


i) 


- X V 


+ goA^- 1 ^ 


fpCLOACVA^V 

2 ° 


1 

+ 




- Pi 1 ’ 0 ) ( A( 2 ) 4 _ 1} + + AW,),) , 


?( 2 ) 


(B. 8 ) 


while the fluid variables get corrections: 


A( n+1 V = ( 1 A ^ oA 

v y n 


(CO 


(h°) 5a Ai 


P 


4 g AM ) » j (^ ) , ( + P)9 , + ^ 2) ) 




AW+ 1 ) 1 ? a = tf 0 E oAr ---^ - 2A oA A^v k A^v k 


d'&or 


- E oA 2 r? 0 AW-i)^ U 2) + 


AW) 


? (D 


+ E oA 2 i? 0 A ( 2 )u; (q 0 A (n X V - C(n-i)) • 

(B.9) 


From here onwards in principle the way would be to solve eqn. (B. 8 ) and find constraints for 
transport coefficients appearing in eqn. (B.7). To solve we would need to plug in the fluid 
variable corrections to all lower orders, along with corrections to lower order constitutive 
relations due to fluid variable corrections. The terms which were zero at equilibrium at lower 
orders will also start to contribute by gaining the fluid variable corrections. Leaving aside 
terms specifically for n = 2, still we would have to deal with a large mess in eqn. (B. 8 ) which is 
analytically not quite tractable. So we leave these expressions at this point for reference. 

Readers are advised that expressions eqns. (B. 8 ) and (B.9) does not contain contributions 
from gravitational and mixed anomaly, and conserved Chern Simons form. Recall that while 
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we set up relations eqn. (3.14), we only used the form of anomalous currents eqn. (2.2) and 

conserved Chern-Simons form eqn. (3.9) to subleading derivative order. At subsubleading 

order, they will receive further gravitational corrections. 
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